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0. Introduction 

Monodromy preserving deformation equations (MPDE) of linear differential operators with rational 
coefficients are known since the beginning of the century [Fu, Schl, G]. Particularly, the famous six Painlcve 
equations are known [G] to be of this type. MPDE were included in the framework of the general theory 
of integrable systems much later, at the end of 70s [ARS, FN1, JMU]; see also [IN]). Many authors were 
inspired by the parallelism between the technique of soliton theory based on the spectral transform and that 
of the MPDE theory based on the monodromy transform. Another issue of this parallelism between soliton 
equations and MPDE is that, in both cases, one deals with certain classes of Hamiltonian systems, namely, 
with infinite-dimensional Hamiltonian structures of evolutionary equations and of their finite-dimensional 
invariant submanifolds in soliton theory, and with remarkable finite-dimensional timc-dcpcndcnt Hamiltonian 
systems in the MPDE theory. 

Recall that one of the first steps in soliton theory was understanding of the Hamiltonian nature of 
the spectral transform as the transformation of the Hamiltonian system to the action-angle variables [ZF] . 
Further development of these ideas was very important for development of the Hamiltonian approach to the 
theory of solitons [FT] and for the creation of a quantum version of this theory. 

In the general theory of MPDE it remains essentially an open question to understand the Hamiltonian 
nature of the monodromy transform, i.e., of the map associating the monodromy data to the linear differential 
operator with rational coefficients. This question was formulated in [FN2] and solved in an example of a 
MPDE of a particular second order linear differential operator. However, the general algebraic properties 
of the arising class of Poisson brackets on the spaces of monodromy data remained unclear. The technique 
of [FN2] seems not to work for more general case. The authors of the papers [AM, FR, KS, Hi] consider 
the important case of MPDE of Fuchsian systems in a more general setting of symplectic structures on 
the moduli space of flat connections (see, e.g., [A]) not writing, however, the Poisson bracket on the space 
of monodromy data in a closed form. MPDE of non-Fuchsian operators and Poisson structure on their 
monodromy data were not considered in these papers. 

In the present paper we solve the problem of computing the Poisson structure of MPDE in the 
monodromy data coordinates for one particular example of the operators with one regular and one irregular 
singularity of Poincare rank 1 

az z 

where U is a diagonal matrix with pairwise distinct entries and V is a skewsymmetric nxn matrix. Recently 
MPDE of this operators proved to play a fundamental role in the theory of Frobenius manifolds [D, Dl]. 
The Poisson structure of MPDE for the operator A coincides with the standard linear Poisson bracket on 
the Lie algebra so(n) 3 V. The most important part of the monodromy data is the Stokes matrix (see the 
definition below). This is an upper triangular matrix S = S(V, U) with all diagonal entries being equal to 1. 
Generically S determines other parts of the monodromy data. It turns out that, although the monodromy 
map 

is given by complicated transcendental functions, the Poisson bracket on the space of Stokes matrices is 
given by very simple degree two polynomials (see formula (3.2) below). The resulting Poisson bracket does 
not depend on U since is involved in the Hamiltonian description of the isomonodromy deformations of the 
operator A(z). The technique of [KS] was important in the derivation of this main result of the present 
paper. 

We hope that this interesting new class of polynomial Poisson brackets and their quantization (cf. [R, 
Ha2]) deserves a further investigation that we are going to continue in subsequent publications. 

The paper is organized as follows: after recalling some basic notations, in section 1.1 we describe the 
monodromy of the operator A(z) around the two singular points; in section 1.2 we present the MPDE for 
this operator. In section 2.1 and 2.2 we describe the related Fuchsian system and its MPDE; in section 2.3 
the Poisson structure on the space of monodromy data of the Fuchsian system is described. In section 3 we 
give the relation between the monodromy data of the two systems and we explicitly calculate the Poisson 
bracket on the space of the Stokes matrices. 



0.1 Basic notations 

Let us consider in the complex domain a differential equation with rational coefficients 

dy 

dz 



^-=A(z)y{z) (0.1) 



where 

v=\"' J \- u-.) 




an Oi2 .. .. din 
«21 0,2n 

* u n l . . .. . . CL nn 

An arbitrary solution y{z) of (0.1) is locally holomorphic but globally multivalued; the poles of A(z) 
are singularities of the solution. Fixing a basis y^ 1 ' , . . . , y^ n ' in the n-dimensional space of solutions we 
construct the fundamental n x n matrix 

^) = (l/ (1) J/ (n) ) 

satisfying the matrix version of (0.1) 

*?M=A(z)Y(z). (0.2) 

dz 

1. Systems with irregular singularity 
1.1 Stokes phenomenon 

In this paper we will concentrate our attention on the linear systems 

^ = (U + -)Y, zeffi, (l.i) 

dz z 

where U is a diagonal nx n matrix with distinct entries m,U2, ■ ■ ■ , u n and V = (vij) € so(n,(D), with 
nonrcsonant eigenvalues Qui, /i2, • • ■ , /Un) (i-e. fj-i — fij (£ Z \ 0). The solutions of the system (1.1) have two 
singular points, and oo. 

• Near the point z = a, fundamental matrix of solutions Yq(z) exists such that 

Y (z) - W{z)z e = [Wo + W lZ , . . .]z e , (1.2) 

where 6 is the diagonalization of V, 6 = W q ~ 1 VWq = diag(/ii, ji2, ■ ■ ■ ,/i„), and VF(z) converges for small 
|z|. Such kind of singularities is called Fuchsian. 

If one continues Y$(z) along a path encircling the point z = 0, the columns of the resulting matrix are 
linear combinations of the columns of Yo(z); hence there exists a matrix Mq such that 

Y Q (z) h-» y (*)M . 

The matrix M is called monodromy matrix around zero. In our case M = exp(2iri8) . 

• At oo the solution has an irregular singularity of Poincare rank 1. This means that it is possible to construct 
a formal series 

T(z) = l+ — + ^ + ... 

z z z 

where V = [Ti, U]+ diagonal, i.e. Ti = (7^) = ( j? ) for i 7^ j, and to define certain sectors &i in which 
a fundamental matrix of solutions Yi exists with asymptotic behavior 

Y i ~Y 00 = T(z)e zU , (1.3) 
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for \z\ — ► oo in Sj. This means that T(z) is the asymptotic expansion of Yie~ zU . 

In different sectors one has different solutions, and this fact is known as Stokes phenomenon. The 
matrices connecting the solutions in different sectors are called Stokes matrices. 

A complete and detailed description of the phenomenon can be found in [BJLI], [Si], [IN],[U]; here we will 
concentrate our attention on the particular operator A(z) = 4- — U — — (see also [D]). 

Following [D] we define an admissible line for the system (1.1) as a line I through the origin on the z-plane 
such that 

Rez(u z - u 3 )\ ze i ^ Vi^j. 

We denote the half-lines 

l + = z : argz = %\) l_ = z : argz = ip — 7T, 

where ip is a fixed real value of the argument. 

The line I lies in the intersection &+ U <5_ of the two sectors 

Si? : xp — ir — e < arg z < ^ + £ 

and 

6l : i\) — e < arg z < V + ^ + e- 

Here e is a sufficiently small positive number. 

Theorem 1.1 : There exists a unique solution Yl(z) analytic in the sector &l with the asymptotic behavior 

yb(«)~Koo; 

the same holds for Yr(z) in &l. 

Proof: See [BJLI]. 

S+ and S- are the Stokes matrices connecting the two solutions in © + , resp. in S_, i.e. 

Y L (z) = Kji(«)S +) z e 6+ 
and 

From the skew-symmetry V T = — V it follows 

o — — *-^+ ■ 

Moreover, one can prove that, given an admissible line, it is possible to order the elements m, i.e., to 
perform a conjugation 

A(z) h^ p- 1 K(z)P, 

where P is the matrix of the permutation in such a way that the Stokes matrix S = 5+ is upper triangular. 

Remark: The full set of monodromy data for the operator A consists of the Stokes matrix S but also of 
the monodromy matrix at the point and of the matrix C connecting the solution (1.2) near zero with a 
solution near the infinity: 

Y (z)=Y L (z)C. 

The monodromy data {S, Mo, C} satisfy certain constraints described in [Dl]. Particularly, 

C- l S T ~ 1 SC = M Q . 



So, in the generic case (i.e., the diagonalizable and nonresonant one) under consideration the diagonal entries 
of Mq e 27rVl , . . . , e 27rv " are the eigenvalues of S T S and C is the diagonalizing transformation for this 
matrix. The ambiguity in the choice of the diagonalizing transformation does not affect the operator A. So, 
the " 2 ~ entries of the Stokes matrix S = S(V; U) can serve as local coordinates near a generic point of 
the space of monodromy data of the operator A (see details in [D], [Dl]). 



1.2 Monodromy Preserving Deformation Equations 

MPDE describe how should the matrix V be deformed, as a function of the "coordinates" Ui , in order to 
preserve the monodromy data. MPDE are the analogue of the isospectral equations in soliton theory. The 
MPDE for the operator A(z) = 4- — U — ¥- are obtained (see [U], [D]) as compatibility equations of the 
system (1.1) with the system 

QY 

i^r {zEi - Vi)Y > 

where Vi — [Pi, Ei\ — adj^.ad^V and (Ei)% = 5f5l. These equations admit the Lax form 

dV 

ss-iw (L4) 

In the generic case (see Remark above) the solution V = V(U) of the MPDE can be locally written in 
implicit form 

S(V;U) = S (1.5) 

for a given constant Stokes matrix S. In other words, the entries of the Stokes matrix can serve as a 
complete system of first integrals of the MPDE (1.4). To explicitly resolve the system (1.5) one has to solve 
an appropriate Riemann-Hilbert boundary value problem. Although this can be explicitly done in a very 
few cases, one can extract certain important information regarding the analytic properties of the solution; 
see more detailed discussions of these properties in [IN] , [JM] , [JMU] , [Si] . 

One can write the MPDE as a Hamiltonian system on the space of the skewsymmetric matrices V 
with the standard linear Poisson bracket for V = (v a b) € so(n): 

{Vab, Vcd} = VadSbc + V bc S ad - V bd S ac - V ac 8 bd . (1.6) 

Indeed, the Lax equation (1.4) can be rewritten as 

dV 



dui 


= {V,Hi(Y,u)}, 


H t = 


_ 1 y- v % 
2 £-?. m- Uj 



for the Hamiltonian function 

1 __ ?; 2 

(1.7) 



In this case, the Poisson bracket is linear but the dynamic of the problem is very complicated; in the following 
we will show how, very much as in the case of isospectral equations, it is possible to find a different coordinate 
system (the entries of the Stokes matrix) in which the dynamic of the evolution is trivial, but the Poisson 
structure is quadratic. The technique developed here consists in building up the monodromy map V — > S 
passing through an auxiliary Fuchsian system. The MPDE for the system (1.1) can be represented also as 
MPDE for an appropriate Fuchsian system 



dx _ \p Aj 

i=\ 



which we shall describe in the next section. The basic idea to construct the Poisson bracket on the space of 
Stokes matrices is to include the map from V G so(n) to S G S into the following commutative diagram of 
Poisson maps 

so(n) — ► S 

I I (1.8) 

A/5 — ► fBl/GL(n,<D) 

where A/% is the space of residues {A{\ of the connection A = ^™=i \- u ^ modulo the action of the 
gauge group %, as we will explain in section 2.2, and dJl/GL(n,(E) is the space of the monodromy data of 
the Fuchsian system (section 2.3), i.e. the space of n-dimensional representations of the free group with n 
generators. 



2. Related Fuchsian system 

2.1 Fuchsian system 

One can relate the system (1.1), with one regular and one irregular singularity to a system with n + 1 
Fuchsian singularities: 



£-e^». <"> 



where 



dX *-^ X — u 

i—l 



B; = -E,[ V+^l 



and 

B 00 = V+h. 

Such a relation is well known in the domain of differential equations, see, e.g. [BJL], [Sch]. 

Now we will briefly describe the monodromy data of the system (2.1). 

In this case Uj is a Fuchsian singular points and, as in (1.2), the general solution near Uj can be expressed 
as 

$ j (X) = W U) (X)(X-u j f\ 

(J) (J) (J) * 1 

where W (A) = W + (X — Uj)W 1 +. . . converges for small |A — Uj\ and Bj = ~^Ej is the diagonalization 

We denote Mj the monodromy matrix along the path "fj encircling the point Uj w.r.t. the basis $oo we 
define in (2.2) below. The matrix Mj is conjugated with the matrix exp(2niBj). 
Also the point oo is Fuchsian; the general solution can be expressed as 

$oo(A)^ <00, (A)(i) Bo \ (2 .2) 

where W °° (A) = W H j \- . . . converges at |A| — > oo and B^ = diag(^ + Hi, ■ ■ ■ , \ + fJ, n ) is the 

diagonalization of B^. Indeed, the following relation holds in the space of the residues : 



n 1 

Y,B i = B 00 = -l + V. 



In this basis the monodromy matrix M^ = —e 2mS . We assume that the loops 71, . . . , j n and 700 are chosen 
in such a way that 

M 1 M 2 ...M n M 00 = l. (2.3) 

2.2 Monodromy Preserving Deformation equations 

We now want to deduce the MPDE for the system (2.1). This amounts to find how can the matrix Bj be 
deformed as function of ui, U2, ■ ■ ■ , u n in order to preserve the monodromy matrices Mi, . . . , M n , M^. The 
answer is given by 

Theorem 2.1 (Schlesinger): // the fundamental solution near infinity is normalized as in (2.2) and Aoo 
is a constant diagonal matrix with nonresonant elements, then the dependence of the Aj on the position of 
the poles of the Fuchsian system 

dX ^ X - u. 
is given, in order to preserve the monodromy, by 



duj ~~ ' i—ii^l Ui—Uj ' 



Proof: it can be found in [Si]. 



Note that system (2.1) does not satisfy the hypotheses of the Schlesinger theorem, because i?oc = 
(V + 5I) is not diagonal. 

In order to apply the Schlesinger theorem it is sufficient to perform the gauge transformation 

B t h+ A t = W^BiWo, (2.4) 

where Wq is the matrix of eigenvectors of V normalized in such a way that 

dWn 

^P- = adE.adu'V. (2.5) 

oui 

Indeed, substituting $ = Wqx, the system (2.1) transforms into 



" Et^t* m 



dX f—' A — u 



and the Schlesinger system follows from the compatibility of (2.6) with 

d X _ M 













dui 


A- Ui 


(See [D]). 
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system 


can 


be 


rewritten 


in the Hamiltonian form 












dA z 
duj 


= {MM 


with the Hamilt 


onians 








«i = - 


y Tr(AjA k ) 
^-^ Ui — Uk 



w.r.t. the linear Poisson bracket 

{{A t )t (A 3 ) c d } = S l3 (S a d (A l ) c b S c b (A 3 ) a d ). (2.7) 

This corresponds to taking, for every Ui, the residue Ai G g[(n,(C) with the natural Poisson bracket on 
g[(n,(D). The residues relative to different singular points commute. In other words (see [KS] , [FR] , [A] ) 
this corresponds to read the matrices Ai as residues of a flat connection (with values in the Lie algebra 
= g[(n,(C)) on the Riemann surface with n + 1 punctures: 

J\ A- 

(in our case uq = oo). On the space of flat connections modulo gauge transformations it is defined the 
Poisson bracket 

r A a , \ A h , M r ab A (Li) — A (v) 

{A (ji),A i/)} = -f c — ^ ^, 

jj, — V 

where / are the structure constants of q w.r.t. the basis {e a } and 

This Poisson bracket gives (2.6). 

Now we can perform the first step in the construction of the map between V and S, that is we pass from 
so(n) to the space A/ 3, where 

n 

A = {V,A 1 ,...,A n \Y / A l = Q} 

i=0 

is the family of the residues of A(X) and 9 is the gauge group. 

Lemma 2.1 :The map V e so(n) i— ► (V,Ai, . . . ,A n ) e A/9 is a Poisson map. (Cf. [Hal], [Hi]) 

Proof: We must compare the Poisson brackets on the two spaces. In so(n) one has the natural coordinates 
{v a b}, with the Poisson bracket (1.6). The natural coordinates in the quotient space A/9 are the traces of 
the products of the matrices Ai, so that we consider the brackets 

{Tr^^TrOM,)} = {(Ai) a b (A k ) b a , (A^A^} = 

= (^)S(^)5{(40«. (^)cl + (AiKiAOiKAkt, {Ai)l}+ 

+ MiAMUAiK, {A l ) d c } + (AkftMiUAiK, {Aj)° d } (2.8a) 

and 

{Tr(AiV),Tr(AjV)} = {(A^V^ {A 3 f d V c d } = 

= (AiMAMiVZrfy + VWUAiK^AiK}. (2.86) 

On A/ 5 by direct calculation, using the bracket (2.7), one obtains 

{TriAiAk^TriAjAt)} = 5 k i r Tr{A i A j A k - A k A 3 A t ) + ^Tr^A^ - A k A t Ai) + 
+ SuTi(A k A j A i - A.AjAk) + ^Tr^^A, - AiA t A k ) = 

= 2[6 k iTr(AiA j A k ) + 6kjTr(A i A l A k ) + 5 t iTr(A k AjA t ) + 6i j Tr(A k AiA i )).(2.9) 



Indeed, A t = —Ei(V + |l) implies 

TrlAiAjAk) = -Tr(A k AjAi) = v ljVok v kl . (2.10) 

On the other hand, Tr (Ai A j) = — u?-, hence 

{Tr( J 4 iJ 4 fe ),Tr( J 4 : ,-A;)} = 4u ifc Uj7{uj fc , ^7} 

= 4:{SkiVijVjkVki - SkjVikVkivu + SuVikVkjVjt + dijVikVkiVu) (2.11) 

where we have used the bracket (1.6). By means of (2.10) it is easy to check that it coincides with (2.9). 

The same can be done for equation (2.8b). Indeed, using the bracket on the Ai matrices and observing 
that 

Tr(A i A j V)=-Tr(VA j A i ) = ^ t>yVfci«j*, 

k^i^j 

one finds 

{Tr(^y),Tr(^F)} = 2Tr(A i A j V). 

On the other hand Tr(AiV) = — X^fe^i v ki' that gi yes 

{Tr(AiV),Tr(AjV)} = 4y]y]vkiVij{vki,vij} = -4^^^^ 

which coincides with (2.8b). Q.E.D 

Lemma 2.2 : The MPDE for the system (1.1) and its related Fuchsian system coincide. 

Proof: It follows immediately from Lemma 2.1 by a straightforward calculation using (2.5), that MPDE 
for the Fuchsian system (2.6) after the gauge transformation (2.4) coincide with (1.4). Actually, one can see 
that the pull back of the Hamiltonian 

^ = y MA 3 Ak) = y TrjBjBk) 

3 ^-^ !1, — Uh U, — Uk 

kjtj 3 k kjtj 3 

is exactely equal to Hj 7 as defined in (1.7). 

2.3 Poisson structure on monodromy data 

In this section we will perform the second step of our construction, that is we will map the Poisson 
structure of A/Q into the space of monodromy data of the Fuchsian system; this is shown in the following 
well-known (see, e.g., [Hi]) 

Theorem 2.2: The monodromy map 

A/5 -» m/SL(n,(E) 

where Wl = {M , Mi, ... , M n \MiM 2 . . ■ M n M = l}, is a Poisson map. 

To actually compute the Poisson bracket on the space of monodromy data, i.e., on the space of n-dimcnsional 
representations of the free group with n generators we will use, following [KS] (Th. 4.2), the following 
technique. We construct the skewsymmetric bracket 



(M t ) a b , {Mj)° d = ztt (MjM^ + (M t M,r d d c b - (M0g(Mj)2 - {Mj) c b {Mi) a d * < j (2.12a) 



(M t )l (M t ) d j = iTT^Mfm (Mf)ffi) . (2.126) 

on the space 9JT of the monodromy matrices. As it was proved in [KS], when restricted to the space 
of representations 9Jt/iSX(n,(D), this bracket defines a Poisson structure on the quotient induced by the 
monodromy map. Observe that the eigenvalues of the matrices Mi are the Casimirs of the Poisson bracket, 
i.e., the functions Poisson commuting with all others (see [KS]). 



3. Poisson structure on the Stokes matrices 



3.1 Connecting the monodromy data of the two systems 

In the previous section we have seen that the space of monodromy data of a Fuchsian system carries a natural 
Poisson stucture. In this section we will show that this structure induces a Poisson bracket on the space of 
Stokes matrices of the related system we studied in chapter 1. To this end we consider the relation between 
the monodromy matrices M\ 1 Mi, . . . , M n of the Fuchsian system and the Stokes matrix S. 
In section 2.1 we claimed that the two systems 

dz z 

and 

d ± = » _^_ 

dX ^ A - u, 
are related, in the sense that, (see Lemma 2.3), the MPDE for the operator A(z) = 4- — U — -j can be 



represented also as MPDE for the operator ^4(A) = -jk- — Ym=i ~^ ± 



dz 

dX /-~ii=l A— u ' 



For a detailed analysis of the transform connecting the two system see [Dl]; here we will concentrate 
our attention on the relation between the monodromy data of the two systems. 

Following Theorem 2.2, we are interested in the quotient of the space of the monodromy data of the Fuchsian 
system w.r.t. the GL(n,(D) conjugations. So, we can choose a particular basis of solution of the system and 
work with the corresponding monodromy matrices. 

Theorem 3.1 : Suppose that yS + S T j is nondegenerate; then there exists a unique basis of solutions (which 

depends on the particular choice of the branchcuts in the complex \~plane) {$ (A)} of the Fuchsian system 
(2.1), such that 

- Near Ui the solution has the behavior 

-m 1 s 



VUi - A a 

- the monodromy matrices are reflections, i.e., going around the singularity Ui the solutions transform 

as 

,«) ,W) „ ,(0 

$ -> $ - 2gij$ 

where G = (gij) = ^ (S + S T ) is the Gram matrix of the following invariant bilinear form w.r.t. the chosen 
basis 

gij = ( * ,$ J :=$ lu-\)& . 

Invariance means that g^ does not depend on A neither on u\, . . . ,u n . 

Proof: Sec [Dl], Th.5.3. 

Remark: $ and Yl are related by the Laplace transform 

2 V^ J -Hi) 
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where 7(j) is a fixed path in the A-plane; analogously for Yr. 



In the {$ (A)} basis the i-th monodromy matrix Mi has the form 



/ 1 




M, 



-2g u -2g 2 



\ 




-2ffr, 



1 / 



\ ... 

2g i0 = 2^ = s^ for i < j. 
This is a reflection w.r.t. the hyperplane normal to the vector <I> . 

The Coxeter identity ([B]) gives 

M 1 M 2 ...M n = -S^S 7 . 

Lemma 3.1: The following relations hold (all the indices are pairwise distinct) 

TriMiMj) =n-A + s 2 i0 
TiiMkMiMjMi) = n - 4 + {s kj - s l3 s lk ) 2 
Tr (M.MjMiMk) = n - 8 + s% + s 2 k + s 2 d + s 2 k + s 2 t + s 2 kl - s i]Sik s jk - 



(3.1a) 
(3.16) 

(3.1c) 



Proof: The fact that the Mi are reflections and that S + S T = 2G geometrically reads into 

where Qjj is the angle between the two hyperplanes normal to $ and $ 

On the other hand, the products MiMj are rotations by the angle 2ctij and this provides the relation 
(3.1a), indeed 



Tr(MjMj-) = n - 2 + 2 cos(2a y ) = n-4 + s 2 j . 



To obtain relation (3.1b) we observe that the product MiMjMi is still a reflection, w.r.t. the mirror normal 
to the vector Mi($ ). This means that the product M k MiMjMi is a rotation by the angle 2/3, where 

-2 cos/3 = [ Mi($ ), $ J = ($ - s y $ , $ ) = s k j - SijSik 

so that T^MkMiMjMi) = n-2 + 2cos(2/3) = n - 4 + (s fej - s l3 s lk ) 2 . 

Finally, (3.1c) can be obtained directly in the case of the 4x4 reflection matrices Mi. Indeed, for 
ordered indices i,j,k,l, the Coxeter identity gives 



where 



M i M j M k M l = -S-^Slu, 



i Sij Si k S%1 

"" I o o 1 Sfe , 

,0001 
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An easy calculation gives the result. 

The same result holds also in dimension n > 4. Indeed, one can observe that, for every n, the product of four 

matrices acts non trivially in the 4-dimcnsional subspacc spanned by the vectors normal to the mirrors of 

the reflections Mi, Mj~, Mi, Mj. It is equal to the identity in the orthogonal complement to the 4-dimensional 

subspace. 

Q.E.D. 

Combining all the above facts, we can conclude our construction proving the following main 

Theorem 3.2: 1. The following formulae 

ITT 

{sik, su} = — (2ski - s ik su) i < k < I (3.2a) 

{sik,Sjk} = y(2sy - s ik Sjk) i < j < k (3.26) 

{sik,Ski} = —{sikSki -2su) i<k<l (3.2c) 

{sik,Sji} = i <k<j <l (3.2d) 

{sik, Sji} = i < j < I < k (3.2e) 

{sife, Sji} = in(sijSki - suSkj) i < j < k <l. (3.2/) 

define a Poisson bracket on the space S of Stokes matrices. 

2. The monodromy map 

so(n) — ► S 

associating the Stokes matrix S G S to the operator A = -4 — U — — , V E so(n), is a Poisson map. 

3. The eigenvalues of S~ 1 S T are the Casimir functions of the Poisson bracket. 

4- The Poisson bracket (3.2) is invariant w.r.t. the action of the braid group B n on the space of braid 
matrices. 

Observe that the Poisson bracket (3.2) does not depend on the times ui, 112, ■ ■ ■ , u n , although the monodromy 
map does. 

Proof: 1. As a first step we explicitly write the restriction of the bracket (2.12) to the space of 
representations. By direct calculation one obtains 

{Tr(M,M fe ),Tr(M J M / )} = {(M,)g(M fe )^ (M^(M^} = 

= {MiY^M^UMkt, (M,)2> + (M,)£(M^{(M fe t (M^}+ 

+ (M fc )£(M,-)S{W)6. (M,)2> + {MktiMttHMit, {Mi)%}. (3.3) 

where we mean summation over repeated indices; using (3.1a), one can rewrite the left hand sides of (3.3) as 

{Tr(M,M fc ), Tr(M,M ; )} = {n - 4 + s 2 lk , n - 4 + s^} = 4« jfc * y {« ifc) a tf }. (3.4) 

Now one has to distinguish between three essentially different cases, in correspondence with the different 
order of the indices. 



• i < k < j < I or i < j < I < k: 

For i < k < j < I all the addenda in the right hand side of (3.3) involve a Poisson bracket of the form 
(2.12) with correctly ordered indices. Here we write explicitly only the first one: 

inTilMiMjMiMk + M l M k M i M j - M l M i M j M k - M l M k M j Mi 
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The others have a similar form, and it is easy to see that they cancel pairwise (the first with the second and 
the third with the fourth). 

The same happens when i < j < I < k, since the only difference is a change of sign in the two last 
elements. Hence it follows 

{Tr(MiM k ), Tr(MjMi)} = i<k<j<l (3.5a) 

{Tr(MiM k ), Tr(MjMi)} = i<j<l<k (3.56) 

Using (3.1b) one immediately obtains equations (3.2d/e) 



• i < j < k < I 

Here the different order of the indices induces a change of sign in the second addendum, which becomes 
equal to the first. Equation (3.3) gives 

{Tr(MiM fe ), Tr(MjMi)} =2iirTr(M i M j MiM k + M l M k MiM J - M l M L M j M k - MiM k MjMi 

=^iTT3 ik Sj l {s i] S k l - SilS k j), 

where the last equality follows from Lemma 3.1. Using eq.(3.4) we obtain immediately eq. (3.2f) 

• i = j < k < I or i < j = k < I or i < j < k = I 

If two indices coincide, for instance i — j < k < I, the other two cases are analogous, wc find 

{Tr(MiM k ),Tr(MiMi)} = {(M 4 )£(M fc )^, (M^M^} = 

= (M^(M^{(M fc )*, (M,)£} + (MOftMOf {(M fc )t (M^}+ 
+ (M fc )*(M i )S{(M i )J, (MM} + (MktiiMOiKMM, {Mi)%}. 

The first and the third addendum cancel, the last is zero (because Mf — 1), and it remains: 



{Tr(M i M fc ),Tr(M i M0} = 2nrUTv{MfM l M k ) - Ty{M % M*M % Mi, 

= 2nr[{n - 4 + s 2 kl ) _ ( n - 4 + s 2 kl + s 2 k s% - 2s k is ik su)} 
= 2iirsi k su(2s k i - SikSu), 

where the second equality follows from (3.1a) and (3.1b). Using (3.4) this leads to (3.2a/b/c). 

2. It follows from the commutativity of the diagram (1.7), where all the arrows are Poisson maps 

3. As we have said above, the eigenvalues of the monodromy matrices are the Casimir functions 
for this Poisson structure. Particularly, applying to M^ we obtain, due to (2.3), the needed statement. 
Practically it is more convenient to use the coefficients of the characteristic polynomial det(S~ 1 S T — /xl) as 
the basic Casimirs. 

4. Recall [D], that the natural action of the braid group B n with n strands on the space of Stokes 
matrices is generated by the following transformations corresponding to the standard generators a\ , . . . , 
c n -i 

Gi : S H^ KiS Ki 
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where the matrix Ki = Ki{S) has the form 

Kjj = l, j = l,...,n; j^i,i + l 

Ka = — Su+i, Ka + i = Ki + u = 1, Ki + u + i = 0. 

Other matrix entries of Ki vanish. According to [D] this action describes the structure of analytic 
continuation of the solutions of MPDE. Our Poisson bracket is obviously invariant w.r.t. analytic 
continuation. 

Q.E.D. 

Example 1. n = 3. In this case the space of Stokes matrices has dimension 3. Denoting x = S12, 
y = S13, z = s 2 3 we obtain, 

, , in , . 

{x,y} = — \2z-xy) 

r 1 in . . 

\ z -, x \ = y( 2 ^~ zx >- 

Our Poisson bracket coincides, within the constant factor — ^, with that of [D]. 

Example 2. n = 4. For convenience of the reader we write here down, omitting the constant factor 
^Y , the Poisson bracket on the six-dimensional space of the Stokes matrices of the form 

(1 p q r' 

1 x y 

1 z 

1, 

{P, q} = (2a; - pq) {x, y} = (2z - xy) 

{P, r} = (2y - pr) {y, z} = (2x - yz) 

U,r} = (2z-qr) {z, x} = (2y - zx) 

{x,p} = (2q-xp) {q, x} = (2p - qx) {r,x} = Q 

{y,p} = (2r- yp) {q,y} = 2(pz - rx) {r, y} = (2p - ry) 

{p, z} = {2, g} = (2r — zq) {r, 2;} = (2q — rz) (3-6) 

The Casimirs of this Poisson bracket are 

C\ = —4 + p + q + r + x + y + z — pqx — pry — qrz — xyz + prxz 

and 

C2 = 6 — 2(+p +q +r +x +y + z ) + 2(—pqx — pry — qrz — xyz) — 2(pqyz + qrxy) +p r +qy +rx 

On the 4-dimensional level surfaces of the Casimirs the Poisson bracket (3.6) induces a symplectic 
structure. These surfaces and the symplectic structures on them are invariant w.r.t. the following action of 
the braid group B4: 

<7i : (p, q, r, x, y, z) ^ (-p, x - pq, y - pr, q, r, z) 

<r 2 : (p, q,r,x, y, z) i-> (q - px,p,r,-x,z - xy, y) 

cr 3 : (p, q, r, x, y, z) ^ (p,r - qz, q,y- xz, x, -z) 
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